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Weyl  points  and  line  nodes  are  three-dimensional  linear  point  and  line  degeneracies  between  two  bands.  In  contrast  to 
two-dimensional  Dirac  points,  which  are  their  lower-dimensional  analogues,  Weyl  points  are  stable  in  momentum  space, 
and  the  associated  surface  states  are  predicted  to  be  topologically  non-trivial.  However,  Weyl  points  are  yet  to  be 
discovered  in  nature.  Here,  we  report  photonic  crystals  based  on  double-gyroid  structures,  exhibiting  frequency-isolated 
Weyl  points  with  complete  phase  diagrams  by  breaking  the  parity  and  time-reversal  symmetries.  Gapless  surface 
dispersions  associated  with  non-zero  Chern  numbers  are  demonstrated.  Line  nodes  are  also  found  in  similar  geometries, 
the  associated  surface  states  forming  flat  dispersion  bands.  Our  results  are  based  on  realistic  ab  initio  calculations  with 
true  predictive  power  and  should  be  readily  realizable  experimentally  from  microwave  to  optical  frequencies. 


Two-dimensional  (2D)  Dirac  points  are  unique  linear- 
dispersion  point  degeneracies  in  2D  periodic  systems.  For 
example,  most  of  the  remarkable  properties  of  graphene  are 
tied  to  the  Dirac  points  at  its  Fermi  level1,2.  In  photonics,  2D 
Dirac  cones  significantly  alter  the  dynamics  of  photons3,4  and 
thus  enable  novel  applications5-7.  Topologically8,  opening  up  the 
2D  Dirac  point  degeneracies  can  generate  topologically  non-trivial 
frequency  gaps  inside  which  unidirectional  ID  edge  states  are  pro¬ 
tected  against  disorder9-12.  Similarly,  it  is  expected  that  if  a  material 
could  be  found  that  exhibits  three-dimensional  (3D)  linear  dis¬ 
persion  relations  between  frequency  (energy)  and  wave  vector  (k), 
it  would  also  display  a  wide  range  of  interesting  phenomena.  Such 
3D  linear  point  degeneracies  between  two  bands  are  called  ‘Weyl 
points’.  Recently,  a  few  proposals13-16  have  been  made  as  to  how 
the  rich  physics  of  Weyl  points17-20  could  be  demonstrated  in  elec¬ 
tronics.  However,  no  observation  of  Weyl  points  has  yet  been 
reported. 

In  this  Article,  we  perform  analytical  and  numerical  studies  to 
predict  a  new  realization  of  Weyl  points  using  3D  photonic  crystals. 
A  key  advantage  of  this  photonic  system,  in  contrast  to  the  proposed 
electronic  systems,  is  the  ease  of  its  eventual  experimental  realization 
and  the  associated  characterization  of  the  Weyl  points.  Specifically, 
we  demonstrate  frequency-isolated  Weyl  points  in  a  perturbed 
double-gyroid  (DG)  photonic  crystal  system.  We  present  complete 
phase  diagrams,  topologically  protected  surface  states,  as  well  as 
frequency- isolated  linear  line  degeneracies  (‘line  nodes’)  and  their 
flat-band  surface  states. 

Before  proceeding,  we  first  point  out  one  intriguing  distinction 
between  the  2D  Dirac  points  and  the  3D  Weyl  points.  2D  Dirac 
cones  are  not  robust;  they  are  only  protected  by  the  product  of 
time-reversal  symmetry  (T)  and  parity  (P)  inversion.  In  two 
dimensions,  the  Dirac  cone  effective  Hamiltonian  takes  the  form 
of  H(k)  =  vxkxax  +  vzkzaz ,  where  v;  are  the  group  velocities  and  cr, 
are  the  Pauli  matrices.  This  form  is  protected  by  PT  (product  of  P 
and  T),  which  requires  H(k)  to  be  real.  Thus,  one  can  open  a  gap  in 
this  dispersion  relation  upon  introducing  a  perturbation,  proportional  to 


that  is  imaginary;  for  example,  even  an  infinitesimal  perturbation  that 
breaks  just  P  or  just  T  will  open  a  gap.  In  contrast,  3D  Weyl  points 
are  topologically  protected  gapless  dispersions  robust  against  any  per¬ 
turbation.  In  three  dimensions,  Weyl  point  dispersions  are  governed 
by  the  Weyl  Hamiltonian  H{ k)  =  vJcx(Jx  +  Vykycry  +  vjizaz.  The  cry 
term  can  exist  only  when  PT  is  broken;  indeed,  this  is  a  necessary  con¬ 
dition  for  the  existence  of  Weyl  points.  Because  all  three  Pauli  matrices 
are  used  in  the  Hamiltonian,  there  is  no  possibility  of  constructing  a 
term  that  can  open  a  gap  in  this  two-band  degeneracy  of  3D  periodic 
systems,  thereby  making  a  single  Weyl  point  absolutely  robust  to  pertur¬ 
bations.  The  only  way  to  eliminate  and  create  Weyl  points  is  through 
pair-annihilations  and  pair-generations  of  Weyl  points  of  opposite  chir¬ 
alities,  which  typically  requires  a  strong  perturbation.  The  chirality 
(c  =  +  l)  of  a  Weyl  point  can  be  defined  as  c  =  sgn(det[v-])  for 
H( k)  =  kfVjjCTj.  It  can  also  be  defined  by  the  Chern  number  of  a  closed 
surface  enclosing  the  single  Weyl  point  in  momentum  space. 

The  organization  of  this  Article  is  as  follows.  Our  starting  point 
for  obtaining  Weyl  points  and  line  nodes  is  a  DG  photonic  crystal 
with  a  band  structure  that  exhibits  a  threefold  degeneracy  (quadratic 
in  all  three  directions)  at  the  Brillouin  zone  centre  (T)  in  an 
otherwise  complete  bandgap.  We  first  obtain  the  line  nodes  by 
applying  a  perturbation  preserving  P  and  T  and  show  the  controlled 
flat  band  surface  states  in  the  surface  Brillouin  zone.  Second, 
we  obtain  the  Weyl  points  by  applying  perturbations  breaking 
either  P  or  T  and  explore  the  phase  diagrams  when  both  P  and  T 
are  broken.  Third,  we  present  the  Chern  numbers  of  the  Weyl 
points  and  the  topologically  non-trivial  surface  states  associated 
with  them.  Finally,  we  discuss  various  opportunities  for 
experimental  realizations. 

Gyroid  photonic  crystals 

The  gyroid,  discovered  by  Schoen21  in  1970,  is  an  infinitely  con¬ 
nected  bi-continuous  triply  periodic  minimal  surface  containing 
no  straight  lines.  Consisting  of  triple  junctions  in  a  body-centred 
cubic  (bcc)  lattice,  the  gyroid  surface  can  be  approximated  by  isosur¬ 
faces  of  g^r)  =  sin(2Trx/fl)cos(2TTy/fl)  +  sin(2'TTy/fl)cos(2'7Tz/a) + 
sin(2Trz/a)cos(2Trx/fl)  (ref.  22),  where  a  is  the  lattice  constant. 
The  space  group  of  a  single  gyroid  (SG)  is  I4j32  (no.  2  1  4)23-25, 
which  lacks  inversion;  adding  its  inversion  counterpart  g(-r) 
gives  the  DG  structure  belonging  to  laid  (no.  230),  which  is  a 
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Figure  1  |  Real-space  unit  cell  and  reciprocal-space  Brillouin  zone  of  the 
gyroid  photonic  crystals,  a,  Real-space  geometry  in  a  bcc  unit  cell  where 
a,  =  (—1, 1, 1)a/2,  a2  =  (1.  —1, 1)a/2  and  a3  =  (1, 1,  — 1)a/2.  The  two 
identical  gyroid  structures  in  red  and  blue  are  high-refractive-index  (n  =  4) 
materials;  they  are  inversion  pairs  with  respect  to  the  origin  (o).  The 
illustrated  air-sphere  of  radius  r  ( r/a=  0.13)  located  at  (1/4,  —1/8, 1/2)a  is 
only  placed  there  when  structural  symmetry  needs  to  be  broken,  b,  The  bulk 
and  (101)  surface  Brillouin  zones  of  the  bcc  lattice.  Weyl  points  and  line 
nodes  investigated  in  this  work  lie  in  the  green  (101)  plane  through  the 
origin  CO  of  the  bulk  Brillouin  zone,  projecting  onto  the  (101)  surface 
Brillouin  zone.  ON  is  along  [101]  and  OH  is  along  [010]  (y).  c,  An 
air-isolated  DG  surface  can  be  formed  by  terminating  the  perturbed  gyroid 
(red)  but  not  the  other  (blue).  The  SG  photonic  crystal  on  top  has  a  large 
complete  bandgap,  as  shown  in  Fig.  2a. 

direct-product  group  of  F4j32  and  inversion.  The  red  gyroid  in 
Fig.  la  is  defined  by  filling  the  inner  space  of  the  isosurface 
{g{ r)  >  1.1)  with  high-refractive-index  material  (^/s  =  n  =  4)  and 
air  otherwise.  The  magnetic  permeability  /z  is  unity  everywhere. 
(These  values  for  the  dielectrics  correspond  to  germanium  and  air 
at  optical  frequencies.)  The  blue  gyroid  is  the  inversion  counterpart 
of  the  red  gyroid  with  respect  to  the  origin;  the  two  gyroids  do  not 
overlap  in  space.  The  band  structures  of  both  the  SG  photonic 
crystal  and  the  DG  photonic  crystal  are  plotted  in  Fig.  2a  in 
orange  and  blue,  respectively.  The  SG  photonic  crystal26  has  a 
32%  complete  bandgap  between  the  second  and  third  bands  from 
0.42  to  0.58  in  normalized  frequencies.  The  DG  photonic  crystal 
band  structure26  contains  a  unique  frequency-isolated  threefold 
degeneracy  among  the  third,  fourth  and  fifth  bands  at  the  centre 
of  the  Brillouin  zone  (T),  which  is  highlighted  by  green  ellipses  in 
Fig.  2a.  The  first  and  second  bands  are  almost  degenerate,  as  are 
the  third  and  fourth  bands,  which  concave  downwards  and  touch 
the  fifth  band,  which  concaves  upwards.  The  threefold  degenerate 
point  is  well  isolated  in  frequency  from  other  states  in  the  band 
structure,  making  it  an  ideal  starting  point  for  applying 
symmetry-breaking  perturbations. 
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Line  nodes  and  their  flat  surface  bands 

The  threefold  degeneracy  of  quadratic  dispersions  at  T  can  be  lifted 
by  breaking  the  74j32  space  group  without  breaking  P  or  T  sym¬ 
metries.  This  is  done  by  replacing  a  part  of  the  gyroid  material 
with  two  air-spheres  (one  on  each  gyroid).  The  first  air-sphere  is 
placed  in  the  red  gyroid  as  illustrated  in  Fig.  la,  and  the  other  is 
its  inversion  counterpart  in  the  blue  gyroid  (not  shown  in 
Fig.  la).  This  perturbation  lifts  the  fifth  band  out  of  the  threefold 
degeneracy  with  the  third  and  fourth  bands  at  T,  as  shown  in 
Fig.  2b.  The  fourth  and  fifth  bands  cross  one  another  linearly, 
forming  a  closed  line  degeneracy  around  the  T  point  in  the  (101) 
plane  through  T,  inside  an  otherwise  complete  frequency  gap. 
The  area  enclosed  by  this  nodal  line  can  be  controlled  by  the 
strength  of  the  perturbations  (the  radii  of  the  air-spheres). 

Similar  to  the  line-node  semimetals27,  the  surface  states  associ¬ 
ated  with  this  line-node  bulk  bandstructure  contain  flat  dispersion 
bands.  We  constructed  an  interface  between  the  DG  and  SG  photo¬ 
nic  crystals  by  removing  only  the  perturbed  gyroid  (Fig.  lc).  The 
suface  states  are  trapped  by  the  pseudo-gap  of  the  DG  and  the  full 
gap  of  the  SG.  We  define  a  termination  parameter  r  (0  <  r<  1) 
to  indicate  the  periodically  equivalent  termination  positions  along 
the  [101]  direction.  r=0  is  set  at  the  origin  of  the  unit  cell,  as 
shown  in  Fig.  lc. 

Figure  3a  shows  one  surface  band  in  the  bulk  pseudo-gap  of  the 
original  DG  photonic  crystal  of  quadratic  point  degeneracy.  When  r 
increases  periodically,  the  surface  dispersion,  at  every  surface  k 
point,  moves  from  the  air  band  (conduction  band)  through  the 
pseudo-gap  to  the  dielectric  band  (valence  band)28.  At  F,  the 
surface  dispersion  is  pinned  into  the  bulk  states  at  the  degeneracy 
point.  The  surface  band  of  r=  0.0,  except  for  the  FI— P  region,  is 
very  flat.  The  high  density  of  states  associated  with  the  flat  surface 
dispersion  is  potentially  useful  for  enhancing  the  light-matter  inter¬ 
actions  at  the  surface.  Even  more  interesting  surface  states  are  shown 
in  Fig.  3b  for  the  line-node  photonic  crystal.  The  nodal  line  bulk 
states  project  onto  the  (101)  surface  Brillouin  zone  as  a  closed  line 
that  separates  the  surface  Brillouin  zone  into  two  disconnected 
areas.  So,  the  surface  dispersions  can  be  flat  bands  in  either  of  the 
two  regions  in  the  Brillouin  zone.  The  green  dispersion  in  Fig.  3b 
has  all  its  frequencies  nearly  degenerate  inside  the  line-node  area, 
while  the  red  dispersion  is  relatively  flat  in  the  rest  of  the 
Brillouin  zone.  The  general  features  of  the  flat  surface  dispersions 
do  not  change  when  the  line-node  photonic  crystal  is  terminated 
by  other  means.  For  example,  when  the  SG  photonic  crystal  is 
replaced  by  air,  one  could  selectively  enhance,  by  changing  the 
surface  terminations,  the  light  emission  of  surface  sources  into 
either  radiative  or  non-radiative  surface  modes  (that  is,  inside  or 
outside  the  light  cone  of  air). 

Phase  diagrams  of  Weyl  points  under  PT-breaking 

In  what  follows,  we  break  the  PT  symmetry  to  obtain  Weyl  points  of 
photons  for  the  first  time.  We  start  by  individually  breaking  P  or  T 
of  the  DG  photonic  crystal,  and  then  consider  the  general  case 
where  P  and  T  are  broken  simultaneously. 

First,  we  break  P  while  preserving  T.  Because  T  maps  a  Weyl 
point  at  k  to  —  k  with  the  same  chirality  (as  velocities  and  cr 
change  signs),  there  must  exist  at  least  two  other  Weyl  points, 
both  of  opposite  chirality,  to  neutralize  the  whole  system.  So,  the 
minimal  number  of  Weyl  points  in  this  case  has  to  be  four.  As 
illustrated  in  Fig.  la,  we  break  P  by  placing  only  one  air-sphere 
on  one  of  the  gyroids  (but  not  the  other)  at  the  middle  point  of 
two  neighbouring  triple  junctions.  Under  this  pure  P-breaking 
perturbation,  two  pairs  of  Weyl  points,  shown  in  Fig.  2c,  emerge 
along  T-N  and  T-H  directions.  The  fact  that  all  the  Weyl 
points  appear  along  high-symmetry  lines  significantly  simplifies 
the  analysis.  There  are  no  other  states  in  the  vicinity  of  the  Weyl 
points’  frequencies. 
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Figure  2  |  Gapless  photonic  band  structures  of  the  DG  photonic  crystals,  a-d,  Band  structures  without  (a)  and  with  (b-d)  perturbations.  All  the  dispersion 
behaviours  close  to  the  degeneracy  points  in  this  figure  can  be  described  well  by  the  low-energy  theory  model  in  Supplementary  Section  SB.  Because  P  and  T 
are  not  broken  at  the  same  time  in  these  photonic  crystals,  k  and  -k  are  degenerate  in  the  band  structures.  A  few  lowest-value  contours  of  the  frequency 
difference  between  the  fourth  and  fifth  bands  are  shown  (hexagonal  insets)  for  each  band  structure  in  the  (101)  plane.  The  contour  spacing  is  0.004  in 
normalized  frequency,  and  '  + '  and  '  are  used  to  label  the  chiralities  of  the  Weyl  points,  a,  The  original  DG  photonic-crystal  band  structure  has  a  threefold 
degeneracy  at  T  among  the  third,  fourth  and  fifth  bands  in  a  pseudo-gap.  The  SG  photonic  crystal  has  a  huge  frequency  gap  covering  the  pseudo-gap 
frequency  region,  b,  We  place  the  two  air-spheres  on  the  two  gyroids:  one  air-sphere  (r/a  =  0.07)  is  located  at  (1/4,  —1/8, 1/2)a,  the  other  is  its  inversion 
symmetric  counterpart.  Under  this  perturbation,  the  fourth  and  fifth  bands  touch  linearly  in  a  closed  line  around  T  in  the  (101)  plane.  The  linear  crossing 
line  is  highlighted  by  a  green  stripe;  this  structure  does  not  yet  exhibit  any  Weyl  points,  c,  We  apply  a  P-breaking  perturbation  by  placing  one  air-sphere 
(r/a=  0.10)  in  one  of  the  gyroids,  but  not  the  other.  Two  pairs  of  Weyl  points  appear  (highlighted  by  green  circles):  one  pair  appears  along  T-H  and  the 
other  along  T-N.  d,  We  apply  a  T-breaking  perturbation  (P  conserved)  by  applying  a  d.c.  magnetic  field  (dimensionless  |B|  =  0.875)  to  the  DG  photonic 
crystal  without  air-spheres.  Only  one  pair  of  Weyl  points  appears  (highlighed  by  the  green  circle),  along  the  direction  of  the  magnetic  field  (T-N). 


We  next  apply  d.c.  magnetic  fields  B  (along  different  directions) 
to  the  original  DG  photonic  crystal  structure  of  Fig.  2a  to  break  T 
while  preserving  P.  We  assume  the  high-index  gyroid  material  is 
gyroelectric  and  use  a  generic  model29  to  describe  its  magnetic 
response.  When  B  is  along  z,  we  assume  the  permittivity  tensor 
takes  the  form  of 


/  eu(|B|)  is12(|B|)  0\ 

e(|B|)  =  — 2e12(|B|)  fill(|B|)  0  (1) 

\  0  Os/ 

where  det(e(|B|))  =  (e2j(|B|)  —  e22(|B|))e  =  s3;this  constant  deter¬ 
minant  condition  ensures  the  dispersions  as  a  whole  do  not  move 
much  in  frequency  with  the  external  d.c.  B  fields.  The  dimensionless 
effective  magnetic  field  intensity  is  defined  as  |B|  =  e12/e  in  this 
Article.  When  the  B  field  is  along  other  directions,  the  correspond¬ 
ing  e  tensor  can  be  obtained  via  coordinate  transformations.  (Note 
that  T-breaking  can  be  implemented  equally  well  via  /Li  for  gyromag- 
netic  materials10).  Under  this  pure  T-breaking  perturbation,  only  a 
single  pair  of  Weyl  points  emerges  along  the  direction  of  the  mag¬ 
netic  field,  as  shown  in  Fig.  2d.  This  is  the  minimum  number  of 
Weyl  points  that  can  exist  with  inversion  symmetry.  These  two 
Weyl  points  are  frequency- degenerate:  P  maps  a  Weyl  point  at  k 
to  —  k  with  the  opposite  chirality  (because  velocities  change  signs). 

Third,  we  apply  both  P-  and  T-breaking  perturbations  continu¬ 
ously,  at  the  same  time,  to  observe  the  phase  transitions  between  the 
two  (II)  Weyl  points  in  the  pure  T-breaking  phase  and  the  four  (IV) 
Weyl  points  in  the  pure  P-breaking  phase.  Interestingly,  different 
magnetic  field  directions  produce  strikingly  different  phase  dia¬ 
grams.  When  B  is  applied  along  T-H,  only  two  phases  exist:  the 
T-breaking  dominated  phase  (II)  and  the  P-breaking  dominated 
phase  (IV).  The  pure  P-breaking  phase,  shown  in  the  contour  plot 


Fig.  2c,  has  four  Weyl  points:  two  with  positive  chiralities  along 
T-H  and  two  with  negative  chiralities  along  T-N.  Applying  a  mag¬ 
netic  field  along  the  T-H  direction  drives  the  two  negative-chirality 
Weyl  points  towards  the  lower  positive  one.  Increasing  the  magnetic 
field  further  annihilates  two  of  the  Weyl  points  of  opposite  chiral¬ 
ities  and  the  system  enters  the  T-breaking  dominated  phase. 


Figure  3  |  Flat  surface  dispersions  of  two  pseudo-gapped  photonic  crystals. 

a,  A  flat  surface  dispersion  at  the  interface  between  the  SG  and  the  original 
DG  of  a  bulk  quadratic  point  degeneracy  as  in  Fig.  2a.  b,  Two  surface 
dispersions  at  the  interfaces  between  the  SG  and  the  line-node  DG  photonic 
crystal  in  Fig.  2b.  They  are  flat  in  two  different  regions  of  the  surface 
Brillouin  zone  separated  by  the  projected  line-node  bulk  states.  The  two 
surface  Brillouin  zones  on  the  right  are  coloured  in  the  regions  where  the 
surface  states  are  flat. 
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Figure  4  |  Phase  diagram  of  Weyl  points  when  the  magnetic  field  is  applied  along  T-N.  a,c,  Two  surface  dispersion  plots  showing  the  entire  band 
dispersions  a)(k)  of  the  third,  fourth  and  fifth  bands  for  two  particular  cases,  where  k  lies  in  the  (101)  surface  Brillouin  zone  in  Fig.  1b.  b,  Each  phase  in  this 
phase  diagram  is  labelled  by  the  number  (Roman  numerals)  of  Weyl  points.  The  phase  transition  line  is  found  by  scanning  the  B  field  amplitude  in  the 
calculation  for  three  different  P-breaking  strengths.  The  pure  P-breaking  phase  (x-axis)  is  highlighted  in  red;  it  is  realizable  at  optical  frequencies  in  essentially 
lossless  materials,  d,  Six  2D  contour  plots  showing  the  locations  of  the  Weyl  points  at  representative  positions  in  the  phase  diagram  (the  contour  spacing  is 
0.002  in  normalized  frequency).  The  contour  plots  share  the  same  colour  as  the  circles  of  their  corresponding  positions  in  the  phase  diagram  in  b. 


which  has  only  two  Weyl  points  along  T-H.  A  detailed  description 
of  this  phase  diagram  is  shown  in  Supplementary  Fig.  SA. 

An  even  richer  phase  diagram,  shown  in  Fig.  4,  appears  when  we 
switch  the  magnetic  field  to  the  T-N  direction.  The  system  under¬ 
goes  two  phase  transitions  from  the  P-breaking  dominated  phase 
(IV)  to  the  T-breaking  dominated  phase  (II).  The  extra  intermediate 
phase  (VI)  of  six  Weyl  points  is  generated  when  one  of  the  negative- 
chirality  Weyl  points  in  the  pure  P-breaking  phase  (red  contour 
plot)  splits  (orange  contour  plot)  into  three  Weyl  points  (blue 
contour  plot)  under  the  increasing  B  field.  Of  these  three  Weyl 
points,  the  original  Weyl  point  flips  its  chirality  from  negative  to 
positive,  and  the  other  two  negative  Weyl  points  move  away  from 
their  creation  position  towards  the  neighbouring  two  positive 
ones.  The  two  pairs  of  Weyl  points  eventually  meet  (green 
contour  plot)  and  annihilate  simultaneously,  leaving  the  remaining 
two  Weyl  points  (purple  contour  plot)  along  the  direction  of  the 
magnetic  field  (T-N).  When  two  Weyl  points  annihilate,  the  band 
dispersions  are  quadratic  along  the  directions  in  which  they  meet. 

In  general,  P-T  phase  diagrams  depend  on  the  form  of  the 
perturbations:  the  d.c.  magnetic  field  can  be  oriented  along 
arbitrary  directions  and,  in  addition,  there  are  many  different 
ways  to  break  the  inversion  symmetry.  We  also  note  that  there  are 
new  topological  phases  outside  of  the  Weyl-point  phases.  For 
example,  by  further  increasing  the  T-breaking  strength  in  phase  II 
of  Fig.  4,  the  two  Weyl  points  could  annihilate  around  the 
Brillouin  zone  boundary  and  enter  a  new  topologically  non-trivial 
phase  of  a  complete  3D  bandgap  that  resembles  the  3D  quantum 
Hall  effect  in  electronics30. 

Chern  numbers  and  non-trivial  surface  states 

Weyl  points  are  topologically  stable  objects  in  the  3D  Brillouin  zone: 
they  act  as  monopoles  of  Berry  flux  in  momentum  space,  and  are 
therefore  intimately  related  to  the  topological  invariant  known  as 
the  Chern  number31.  The  Chern  number  can  be  defined  for  a 
single  bulk  band  or  a  set  of  bands,  where  the  Chern  numbers  of 
the  individual  bands  are  summed,  on  any  closed  2D  surface  in 
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the  3D  Brillouin  zone.  The  difference  of  the  Chern  numbers 
defined  on  two  surfaces,  of  all  bands  below  the  Weyl  point  frequen¬ 
cies,  equals  the  sum  of  the  chiralities  of  the  Weyl  points  enclosed 
between  the  two  surfaces18.  In  Fig.  5a  we  illustrate  the  Chern 
numbers  of  the  2D  planes  perpendicular  to  T-H(y);  these  are 
closed  surfaces  in  the  3D  periodic  Brillouin  zone.  The  first 
contour  plot  (red),  of  a  pure  P-breaking  phase  in  Fig.  4,  is  used  as 
an  example.  The  Chern  number  of  all  lower  bands  on  the  plane  is 
plotted  in  light  blue  on  the  left  side.  The  Chern  number  vanishes 
when  the  plane  is  outside  the  Weyl  points.  It  changes  by  1  when 
the  plane  moves  across  one  Weyl  point  and  changes  by  2  when 
the  plane  crosses  two  Weyl  points  of  the  same  chirality. 

As  a  result  of  the  non-zero  Chern  numbers  of  all  the  lower  bulk 
bands  (first,  second,  third  and  fourth),  there  are  topologically  pro¬ 
tected  gapless  chiral  surface  states  inside  the  bandgap  (between 
fourth  and  fifth  bands)  away  from  the  Weyl  points.  This  is  a 
higher-dimensional  generalization  of  the  protection  of  one-way 
chiral  edge  states  (ID)  by  the  non-zero  Chern  numbers  of  the  2D 
bulk  bands13,32.  An  example  of  the  non-trivial  gapless  surface  dis¬ 
persion,  along  a  line  cut  in  the  2D  surface  Brillouin  zone  in 
Fig.  5a,  is  plotted  in  Fig.  5b.  A  surface  mode  profile  is  also  shown. 
Because  T  is  not  broken  in  this  example,  the  surface  dispersions 
are  degenerate  between  k  and  —  k  while  the  Chern  number  flips 
sign.  Similar  to  the  ‘Fermi  arcs’  in  Weyl  semimetals13,18,  at  the  fre¬ 
quency  of  the  Weyl  points,  the  surface  dispersions  connect  the 
Weyl  points  (of  opposite  chiralities)  as  line  segments  that  do  not 
close  in  the  surface  Brillouin  zone.  The  above  analysis  on  Chern 
numbers  and  surface  modes  applies  for  surfaces  terminated  along 
other  directions  and  for  other  Weyl-point  phases  as  well. 

Proposals  for  experimental  realizations 

The  theoretical  predictions  discussed  above  can  be  readily  realized 
using  many  available  materials  and  fabrication  techniques,  and 
the  bulk  and  surface  dispersions  can  be  verified  by  transmission 
experiments33,34.  The  line  nodes  and  the  Weyl-point  phases  of 
pure  P-breaking  perturbations,  without  T-breaking,  can  be  realized 
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P-broken  DG  photonic  crystal  SG  photonic  crystal 


Figure  5  |  Topological  surface  states  of  a  Weyl-point  photonic  crystal. 

a,  The  non-zero  Chern  number  of  the  plane  moving  along  ky  in  the  3D  bulk 
Brillouin  zone.  The  Chern  number  is  summed  over  the  bulk  bands  below  the 
Weyl-point  frequency,  b,  An  example  of  a  topologically  protected  surface 
state  of  a  Weyl-point  photonic  crystal.  The  projected  bulk  states  are  plotted 
in  blue  in  the  surface  bandstructure,  and  the  gapless  surface  state  dispersion 
relation  is  plotted  in  red.  The  spatial  field  distribution  (|E|2)  of  one  surface 
mode,  at  the  centre  of  the  k-axis,  is  plotted  with  the  contour  outlines  of  the 
refractive  index  profile.  The  surface  state  is  localized  at  the  interface 
between  the  DG  photonic  crystal  and  the  SG  photonic  crystal. 

from  optical  to  microwave  frequencies.  The  Weyl-point  phases 
involving  T-breaking  can  be  realized  with  gyromagnetic  materials 
at  microwave  frequencies10. 

The  gyroid  structures  can  be  made  by  either  bottom-up  or  top- 
down  approaches.  In  nature,  gyroid  geometries  are  found  in  butter¬ 
fly  wings35,  lipid  mesophases36,  liquid  crystals37  and  block  co-poly- 
mers38'39.  At  optical  wavelengths,  germanium  can  either  be  directly 
synthesized40  into  DG  strucutures  or  backfilled41  into  DG  templates 
made  by  self-assembly.  Top-down  3D  fabrications  such  as  interfer¬ 
ence  lithography42,  direct  writing43,  angled-etching  (for  the  SG)44,45, 
membrane-stacking46  and  3D  printing  are  all  good  candidates  for 
making  the  gyroid  structures. 

Conclusions 

In  summary,  this  article  provides  a  detailed  design  of  the  minimum 
number  of  frequency-isolated  Weyl  points  and  line  nodes  along  the 
high-symmetry  lines  and  planes  in  the  Brillouin  zone  of  a  DG 
photonic  crystal  system.  Complete  P-T  phase  diagrams  are  investi¬ 
gated  by  means  of  annihilations  and  creations  of  Weyl  pairs.  The 
topologically  non-trivial  surface  states  of  Weyl-point  photonic  crys¬ 
tals  and  flat  surface  dispersions  of  line-node  photonic  crystals  are 
presented.  This  work  may  open  doors  to  new  paradigms  in  photo¬ 
nics:  topologically  protected  2D  chiral  surface  states  realizable  at 
optical  frequencies,  radiation-controllable  photonic  density  of 
states  enhanced  by  flat  surface  dispersions,  possible  new  3D  topolo¬ 
gical  phases  of  complete  frequency  gaps  by  annihilations  of  the 
Weyl  points,  novel  transmission  properties  and  other  unconven¬ 
tional  phenomena  associated  with  the  density  of  states  and  dis¬ 
persion  relations  close  to  the  degeneracy  points. 
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